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Abstract.  A  genetic  and  morphologic  survey  of  Anopheles  darlingi  populations  collected  from  seven  countries  in 
Central  and  South  America  was  performed  to  clarify  the  taxonomic  status  of  this  major  malaria  vector  species  in  the 
Americas.  Population  genetics  was  based  on  three  techniques  including  isozyme,  random  amplified  polymorphic 
DNA-polymerase  chain  reaction  (RAPD-PCR),  and  internal  transcribed  spacer  2  (ITS2)  markers.  The  results  of  the 
isozyme  analysis  indicated  moderate  differences  in  the  allele  frequencies  of  three  putative  loci  (glutamate  oxalaoace- 
tate  transaminase- 1,  isocitrate  dehydrogenase- 1,  and  phosphoglucomutase)  of  the  31  analyzed.  No  fixed  electromorph- 
ic  differences  separated  the  populations  of  An.  darlingi,  which  showed  little  genetic  divergence  (Nei  distances  = 
0.976-0.995).  Fragments  produced  by  RAPD-PCR  demonstrated  evidence  of  geographic  partitioning  and  showed  that 
all  populations  were  separated  by  small  genetic  distances  as  measured  with  the  1  -  S  distance  matrix.  The  ITS2 
sequences  for  all  samples  were  identical  except  for  four  individuals  from  Belize  that  differed  by  a  three-base  deletion 
(CCC).  The  morphologic  study  demonstrated  that  the  Euclidean  distances  ranged  from  0.02  to  0. 14,  with  the  highest 
value  observed  between  populations  from  Belize  and  Bolivia.  Based  on  these  analyses,  all  the  An.  darlingi  populations 
examined  demonstrated  a  genetic  similarity  that  is  consistent  with  the  existence  of  a  single  species  and  suggest  that 
gene  flow  is  occurring  throughout  the  species’  geographic  range. 


Anopheles  (Nyssorhynchus)  darlingi  Root  1926  is  the 
most  efficient  malaria  vector  in  the  New  World.'  In  South 
America,  its  distribution  stretches  from  Colombia  to  north¬ 
eastern  Argentina.  In  Mexico,  it  is  found  as  far  north  as 
Tabasco  and  Chiapas  States,  and  in  Central  America,  it  oc¬ 
curs  in  Belize,  Guatemala,  Honduras,  and  occasionally  in  El 
Salvador.'^  However,  this  species  has  never  been  officially 
reported  from  Nicaragua,  Costa  Rica,  or  Panama,  resulting 
in  an  apparent  discontinuity  in  its  distribution  between  Cen¬ 
tral  and  South  America." 

Transmission  of  human  malaria  by  An.  darlingi  occurs 
mainly  in  warm  and  humid  forests  along  rivers  in  lowland 
areas  (less  than  550  meters  above  sea  level),  sometimes  in 
coastal  and  piedmont  ecoregions  with  high  rainfall  and  trop¬ 
ical  forests."*  Anopheles  darlingi  is  a  very  serious  malaria 
vector  throughout  its  geographic  range  because  of  its  asso¬ 
ciation  with  humans  and  its  high  rates  of  infection  with  the 
malaria  parasites  Plasmodium  falciparum,  P.  vivax,  and  P. 
malariae.^  It  is  the  most  anthropophilic  and  endophagic 
anopheline  in  the  Americas,  factors  that  obviously  account 
for  much  of  its  potency  as  a  major  malaria  vector."”*  As  a 
result  of  complex  sociologic,  economic,  and  biologic  factors, 
malaria  has  resurged  in  Brazil."*  Anopheles  darlingi  is  the 
primary  vector  of  human  malaria  in  the  Amazon  Basin,"  and 
is  the  major  contributor  to  the  present  resurgence  of  disease 
in  South  and  Central  America. 

Anopheles  darlingi  was  first  described  by  Root"’  in  1926 
from  Rio  de  Janeiro  in  Brazil.  In  1937,  from  samples  col¬ 
lected  near  Sao  Paulo  (Brazil),  Galvao  and  others"  described 
a  variety  of  An.  darlingi  named  paulistensis  based  on  adult 
and  egg  characteristics.  The  latter  characters  seem  to  be  at 
least  quite  variable  in  Brazil'-""  and  Venezuela.'-*  Studies  on 
the  biology  of  An.  darlingi  also  showed  considerable  varia¬ 


tion  in  Brazil,*  Suriname,"**  "'  French  Guiana,""  Guyana,"-" 
Colombia,""*  and  Venezuela.""  In  addition,  variation  has  been 
described  in  its  morphology,"*""  chromosomes,"*  "’  and  iso¬ 
zymes.'*’*'-""  Some  of  these  reports  have  suggested  that  An. 
darlingi  may  be  a  species  complex.'’-"® 

In  1988,  Linthicum"  asserted  that  populations  of  An.  dar¬ 
lingi  collected  throughout  its  geographic  range  showed  little 
interpopulation  variation  among  the  diagnostic  morphologic 
characteristics.  Therefore,  An.  darlingi  was  defined  as  a 
monotypic  subgroup  and  placed  in  the  Argyritarsis  section. 
However,  a  recent  phylogenetic  analysis  (Danoff-Burg  JA, 
Conn  IE,  unpublished  data)  using  morphologic  and  molec¬ 
ular  data  suggested  that  the  Argyritarsis  section  is  paraphy- 
letic  with  respect  to  the  Albimanus  section  and  An.  darlingi 
has  been  placed  in  a  basal  grade  in  the  Nyssorhynchus  sub¬ 
genus  (Figure  1). 

In  spite  of  these  reports,  few  comparative  studies  of  An. 
darlingi  populations  have  been  undertaken.  The  present 
study  was  based  on  populations  of  An.  darlingi  from 
throughout  much  of  its  geographic  range  in  Central  America 
and  South  America  (Figure  2).  In  Brazil,  An.  darlingi  was 
collected  in  different  states  (Table  1),  including  the  type- 
locality  of  the  variant  paulistensis  (from  Sao  Paulo  State). 

The  approach  we  took  to  evaluate  the  taxonomic  status  of 
An.  darlingi  was  based  on  morphometric  and  genetic  analy¬ 
ses  using  isozyme,  random  amplified  polymorphic  DNA- 
polymerase  chain  reaction  (RAPD-PCR),  and  internal  tran¬ 
scribed  spacer  2  (ITS2)  markers.  These  genetic  techniques 
have  proven  to  be  excellent  for  distinguishing  species  of 
insects""-""  that  otherwise  are  very  difficult  to  separate  on 
the  basis  of  morphologic  features.  Our  main  objectives  in¬ 
cluded  1)  analysis  of  the  differentiation  and  variability  of 
An.  darlingi,  2)  comparison  of  An,  darlingi  populations  with 
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Figure  1.  Basal  portion  of  Nyssorhynchus  subgenus  phylogeny 
based  on  data  from  complete  sequences  of  NADH  dehydrogenase 
subunit  2,  NADH  dehydrogenase  subunit  6,  cytochrome  oxidase 
subunit  I,  and  internal  transcribed  spacer  2,  as  well  as  egg,  larval, 
pupal,  and  adult  morphologic  characteristics.^^  Species  in  italics  are 
those  that  were  postulated  to  be  in  the  Albimanus  section,  whereas 
those  in  roman  font  were  said  to  be  in  the  Argyritarsis  section.^’^*  ™ 
Relationships  among  the  terminal  species  in  the  subgenus  were  sim¬ 
plified  to  demonstrate  the  basal  placement  of  Anopheles  darlingi. 


the  variant  paulistensis  described  in  the  literature,  and  3) 
comparison  of  the  morphologic  and  genetic  profiles  of  An. 
darlingi  populations  with  other  Anopheles  species  used  as 
outgroups. 

MATERIALS  AND  METHODS 

Mosquito  populations.  From  1990  to  1996,  we  collected 
samples  of  An.  darlingi  from  seven  countries  including  Be¬ 
lize  in  Central  America,  and  Bolivia,  Brazil,  Colombia, 
French  Guiana,  Peru,  and  Venezuela  in  South  America  (Fig¬ 
ure  2  and  Table  1).  The  countries  were  chosen  to  provide  a 
spatial  representation  of  the  geographic  distribution  of  this 
species  and  to  examine  samples  from  the  type-locality  of  the 
variant  paulistensis.  Populations  (Per  1-5)  from  Pera  (Iqui- 
tos,  3°46'S,  73°15'W)  were  used  only  for  the  RAPD-PCR. 
In  the  morphologic  study,  populations  of  An.  darlingi  from 
northern  Colombia  (Bete,  Departamento  del  Chocd,  5°45'N, 
77°W)  were  also  analyzed. 

Samples  used  were  primarily  single  individuals  from  iso¬ 
female  progeny  broods,  wild-caught  adults,  and  occasionally 
larvae  or  pupae  reared  individually  to  the  adult  stage.  For 
the  latter,  fourth  instar  larval  and/or  pupal  exuviae  were  pre¬ 
served  and  each  specimen  was  recorded,  identified  to  spe¬ 
cies,  and  kept  frozen  at  —  70°C. 

Other  Anopheles  species  were  included  in  the  isozyme, 
RAPD-PCR,  and  morphologic  analyses  as  outgroup  taxa.  We 
used  An.  (Nys.)  argyritarsis  Robineau-Desvoidy  and  An. 
(An.)  pseudopunctipennis  Theobald  for  the  isozyme  analysis. 
An.  (Nys.)  nuneztovari  Gabaldon  for  the  RAPD-PCR,  and 


Figure  2.  Collecting  areas  of  Anopheles  darlingi  in  seven  coun¬ 
tries  throughout  the  geographic  distribution  of  this  species.  Each 
symbol  (squares  and  circles)  corresponds  to  a  collection  of  multiple 
samples.  Dark  squares  are  samples  used  in  genetic  analyses,  white 
squares  in  morphologic  studies,  and  dark  circles  in  both  genetic  and 
morphologic  studies.  ALB  =  Fortuna  de  Albarical  (Amazonas,  Ven¬ 
ezuela);  ARA  =  Araraquara  (Sao  Paulo,  Brazil);  ARI  =  Aripao 
(Bolivar,  Venezuela);  BAH  =  Itabela  (Bahia,  Brazil);  BAR  =  Pereira 
Barreto  (Sao  Paulo,  Brazil);  BEL  and  BHZ  =  Belize  (Central  Amer¬ 
ica);  BEM  =  Belem  (Pard,  Brazil);  BOL  =  Beni  (Bolivie);  CAP  = 
Capanema  (Para,  Brazil);  CAS  =  Caicara  (Bolivar,  Venezuela); 
CMB  =  Costa  Marques  (Rondonia,  Brazil);  COL  =  Bete  (Choco, 
Colombia);  CUI  =  San  Antonio  de  Leverger  (Mato  Grosso,  Brazil); 
DOU  =  Rio  Jacare  Pepira  (Sao  Paulo,  Brazil);  ELI  =  El  Juval 
(Trujillo,  Venezuela);  IQ  =  Iquitos  (Pern);  ITU  =  Ituxi  River  (Ama¬ 
zonas,  Brazil);  MAN  =  Manaus  (Amazonas,  Brazil);  MAR  =  Mar- 
ipasoula  (French  Guiana);  MG  =  Mato  Grosso  (Brazil);  OCA  = 
Ocamo  (Amazonas,  Venezuela);  PEX  =  Peixoto  de  Azevedo  (Mato 
Grosso,  Brazil);  PIN  =  El  Pinal  (Bolivar,  Venezuela);  PR  =  Velha 
Timboteua  (Pard,  Brazil);  SIM  =  San  Juan  de  Manapiare  (Amazon¬ 
as,  Venezuela);  TRP  =  Trou-Poisson  (French  Guiana);  URU  =  Um- 
curi  (Pard,  Brazil);  YUC  =  La  Yuca  (Bolivar,  Venezuela). 


four  species  in  the  Nyssorhynchus  subgenus;  An.  albitarsis 
Lynch  Arribalzaga,  An.  argyritarsis.  An.  braziliensis  (Cha¬ 
gas’),  and  An.  marajoara  (Galvao  and  Damascene)  for  the 
morphologic  analysis.  When  relationships  among  21  other 
Nyssorhynchus  taxa  were  reconstructed  using  only  ITS2,  An. 
darlingi  was  monophyletic  and  included  in  a  basal  unre¬ 
solved  polytomy.  Therefore,  no  outgroups  were  used  in  the 
present  study. 

Isozyme  electrophoresis.  Isozymes  were  separated  by 
horizontal  starch  gel  electrophoresis.  Experimental  proce¬ 
dures  are  detailed  in  the  publication  by  Manguin  and  oth- 
ers."**  A  total  of  35  enzyme  systems  was  screened  on  two 
different  buffers:  the  morpholine  buffer  system"'^  (Morph,  pH 
6.1),  and  the  Tris-citrate  buffer  system"  (TCss,  pH  6.7).  Of 
the  35  enzyme  systems  tested,  23  showed  good  allehc  res¬ 
olution,  including  31  putative  loci  (Table  2).  For  each  locus, 
the  most  frequent  electromorph  was  designated  the  100  al¬ 
lele  and  all  other  alleles  were  measured  relative  to  it.  Elec¬ 
tromorph  genotype  frequencies  were  used  as  input  for  BIO- 
SYS-l."*^  Analysis  of  each  population  included  computation 
of  allele  frequencies,  heterozygosity  per  locus,  additional 


Table  1 

Geographic  information  on  the  sample  sites  of  Anopheles  darlingi 
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Table  2 

Electrophoretically  detected  enzyme  systems  of  Anopheles  darlingi 


Enzyme  system 

E.C.  number’* 

Symbol 

No,  of  locit 

Buffert 

1.  Adenylate  kinase 

2.7.4.3 

AK 

1 

Morph 

2.  Aldehyde  oxidase 

1.2.3.1 

AO 

1 

Morph 

3.  Arginine  kinase 

2.7.3.3 

ARGK 

1 

TCss 

4.  Esterase 

3.1. 1.1 

EST 

2 

Morph 

5.  Fumarase 

4.2.1.2 

FUM 

1 

TCss 

6.  Glvcerol  dehydrogenase 

1.1.1.72 

GCD 

1 

Morph 

7.  Glutamate  oxaloacetate  transaminase 

2.6.1. 1 

GOT 

2 

Morph 

8.  a-Glycerophosphate  dehydrogenase 

1.1. 1.8 

GPDH 

1 

Morph 

9.  Glyceraldehyde-3-phosphate  dehydrogenase 

1.2.1.12 

G3PDH 

1 

Morph 

10.  Glutathione  reductase 

1. 6.4.2 

GR 

2 

TCss 

1 1 .  fl-Hydroxyacid  dehydrogenase 

1.1.1.30 

HAD 

1 

Morph 

12.  Hexokinase 

2.7.1. 1 

HK 

2 

Morph 

13.  Isocitrate  dehydrogenase 

1.1.1.42 

IDH 

2 

TCss 

14.  Leucine  amino  peptidase 

3.4.11.1 

LAP 

2 

TCss 

15.  Malate  dehydrogenase 

1.1.1.37 

MDH 

2 

Morph 

16.  Malic  enzyme 

1.1.1.40 

ME 

1 

Morph 

17.  Mannose-6-phosphate  isomerase 

5.3. 1.8 

MPI 

1 

Morph 

18.  6-Phosphog!uconate  dehydrogenase 

1.1.1.44 

6PGD 

1 

Morph 

19.  Phosphoglucose  isomerase 

5.3. 1.9 

PGI 

1 

TCss 

20.  Phosphoglucomutase 

5. 4.2.2 

PGM 

1 

Morph 

21.  Pyruvate  kinase 

2.7.1.40 

PK 

2 

Morph 

22.  Triose  phosphate  isomerase 

5.3. 1.1 

TPI 

1 

TCss 

23.  Xanthine  dehydrogenase 

1.2.1.37 

XDH 

1 

TCss 

Total 

31 

*  Enzyme  commission  number, 
t  Number  of  scorable  bands  per  phenotype. 

t  Refers  to  electrophoresis  buffer:  Morph  =  morpholine;  TCss  =  Tris-citrate. 


measures  of  genetic  variability,  and  a  test  for  conformance 
to  Hardy- Weinberg  equilibrium  at  single  loci  by  chi-square 
analysis.  Differentiation  among  the  populations  was  mea¬ 
sured  by  F-statistics.  Nei’s'**  unbiased  genetic  distances  and 
identities  and  Rogers'*^  genetic  distances  were  clustered  by 
the  unweighted  pair  group  method  using  arithmetic  average 
(UPGMA)  to  produce  the  phenogram. 

Random  amplified  polymorphic  DNA-polymerase 
chain  reaction.  Experimental  procedures  can  be  found  in 
the  publications  by  Wilkerson  and  others.^’^’'**  Sixty  deca- 
mer  primers  (Operon  sets  A,  B,  C;  Operon  Technologies, 
Alameda,  CA)  were  screened  using  three  individuals  from  a 
single  family  of  An.  darlingi  from  Peixoto  de  Azevedo  (Bra¬ 
zil).  From  these,  a  subset  of  eight  primers  that  produced 
consistent  bands  that  could  be  scored  in  each  of  the  three 
specimens  was  chosen  to  analyze  the  entire  sample  of  124 
individuals  from  13  populations  (Table  1).  Five  of  the  eight 
gave  consistent  results  for  all  specimens  in  the  study.  In 
contrast,  13  individuals  of  An.  nuneztovari  from  two  popu¬ 
lations  were  analyzed  using  the  same  five  primers.  Anoph¬ 
eles  nuneztovari  is  placed  in  a  sister  clade  to  the  one  in¬ 
cluding  An.  darlingi.  The  resultant  data  set  was  formatted  as 
described  by  Black‘d  and  a  1  —  S  distance  matrix  was  gen¬ 
erated  using  the  similarity  option  in  the  RAPDPLOT  (FOR¬ 
TRAN  programs  for  the  analysis  of  RAPD-PCR  markers  in 
populations;  Black  IV,  WC,  Colorado  State  University,  Fort 
Collins,  CO)  program.  The  formula  is  derived  from  the  Nei 
and  LP®  similarity  index:  S  =  2Nab/(Na  +  Ng)  where  N^b 
are  the  fragments  that  two  individuals  share  in  common  and 
Na  and  Ng  are  the  number  of  fragments  in  individuals  A  and 
B,  respectively.  The  matrix  was  analyzed  in  PHYLIP  3.5C 
(Phylogeny  inference  package;  Felsenstein  J,  University  of 
Washington,  Seattle,  WA)  using  the  NEIGHBOR  program 


with  the  UPGMA  option  and  a  phenogram  was  produced 
with  DRAWGRAM,  also  in  PHYLIP  3.5C  RAPDBOOT^' 
was  used  to  generate  1(X)  pseudoreplicate  distance  matrices, 
which  were  collapsed  to  form  100  trees  with  UPGMA.  The 
bootstrap  consensus  tree  was  derived  from  the  100  UPGMA 
trees  with  the  CONSENSUS  program  in  PHYLIP  3.5C. 

Sequencing  of  ITS2.  Total  DNA  was  extracted  from  spec¬ 
imens  using  the  extraction  technique  of  Collins  and  others.’^ 
The  PCR  temperature  profile  consisted  of  37  cycles  of  1  min 
at  95°C,  1  min  at  45°C,  and  1.5  min  at  72°C.  Primers  specific 
for  the  ITS2  region^^  annealed  to  the  conserved  regions  of 
the  5.8S  and  28S  ribosomal  subunits  that  flank  this  region. 
All  13  specimens  were  manually  sequenced  and  represen¬ 
tatives  from  each  of  the  12  sites  were  resequenced  using  an 
ABI  373  automated  sequencer  using  the  ABI  PRISM  Dye 
terminator  cycle  sequencing  kit  (Perkin  Elmer,  Norwalk, 
CT).  Contig  assembly  and  related  sequence  analysis  was  per¬ 
formed  using  Sequencher  3.0  (Gene  Codes  Corp.,  Ann  Ar¬ 
bor,  MI).  The  GenBank  accession  number  of  the  ITS2  for 
An.  darlingi  populations  from  South  America  is  U92337, 
and  the  number  for  populations  from  BeUze  is  AF051274 
(Danoff-Burg  JA,  Conn  JE,  unpublished  data). 

Morphology.  Selected  wing,  leg,  and  head  characteristics 
were  measured  and  analyzed  on  305  adult  females  of  An. 
darlingi  from  17  populations  (Figure  2  and  Table  1)  using  a 
dissecting  microscope.  The  head  characteristic  was  the 
length  of  the  proboscis  and  those  of  the  wings  and  legs  are 
illustrated  in  Figure  3  A  and  B,  respectively.  The  morpho¬ 
logic  terms  and  abbreviations  used  in  this  figure  follow  those 
of  Harbach  and  Knight,^'*-^^  Wilkerson  and  Peyton,^®  and  Pey¬ 
ton  (unpublished  data).  The  Program  Statistical'^  was  used 
for  statistical  analysis  among  populations  with  N  >10.  To 
measure  the  closeness  of  An.  darlingi  populations  based  on 
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FlGl^  3.  Morphologic  wing  and  leg  characteristics  of  Anopheles  darlingi.  A,  wing  types  based  on  patterns  of  dark  and  pale  spots  on  the 
costa.  BP  +  PHP  -  basal  pale  +  prehumeral  pale;  PHD  =  prehumeral  dark;  h  =  humeral  crossvein;  HP  =  humeral  pale-  HD  =  humeral 
dark;  PSP  =  presector  pale;  PSD  =  presector  dark;  SP  =  sector  cale;  SD  =  sector  dark;  PRSD  =  proximal  sector  dark-  aqp  =  a-'-essory 
sector  pale;  DSD  =  distal  sector  dark;  SCP  =  subcostal  pale;  PD=  preapical  dark;  PP  =  preapical  pale;  AD  =  apical  dark;  A^P  =  apical  pale 
(wing  spot);  ii  rj  -  radius  one  and  cell  minus  radial  sector.  B,  legs  showing  the  characteristics  measured.  Fe-I  =  forefemur;  Ta-1  = 
foretarsomere  1;  pale  spot  (PS-1,)  on  Ta-I,;  Ta-Ij  =  foretarsomere  2;  pale  spot  (PS-lj)  on  Ta-Ij;  Ta-Ij  =  foretarsome  3;  pale  spot  (PS-1,)  on  Va- 
I,;  Ta-IIIj  =  hindtarsomere  2;  pale  spot  (PS-lIIj)  on  Ta-Illj. 


linear  continuous  characteristics  such  as  morphologic  ones, 
the  Euclidean  distances  among  populations  were  estimated 
by  cluster  analysis.^® 

RESULTS 

Isozyme  analysis.  An  isozyme  comparison  of  31  loci  (Ta¬ 
ble  2)  among  the  populations  of  An.  darlingi  indicated  that 
mean  heterozygosity  ranged  from  0.063  to  0.122  (Table  3) 
with  a  mean  ±  SE  of  0.1  ±  0.023  across  all  populations. 

The  F-statistic  (Fst),  a  measure  of  the  amount  of  differ¬ 


entiation  among  subpopulations,  showed  a  low  mean  value 
of  0.102  and  the  mean  index  of  fixation  of  individuals  rel¬ 
ative  to  the  total  of  subpopulations  (Fjg)  had  a  moderate 
mean  value  of  0.039.  Three  Fs^  values  with  moderate  dif¬ 
ferentiation,  were  found  for  glutamate  oxaloacetate  trans- 
atninase-1  (Got-1)  (0.170),  isocitrate  dehydrogenase- 1  {Idh- 
1)  (0.170),  and  phosphoglucomutase  (Pgni)  (0.204). 

There  were  five  different  alleles  for  Got-1.  Populations  of 
An.  darlingi  from  Sao  Paulo  (Brazil)  had  a  62%  frequency 
for  allele  100  and  37%  for  allele  85,  whereas  all  the  other 
populations  had  a  frequency  more  than  80%  for  allele  100 
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Table  3 

Measures  of  genetic  variatoin  of  Anopheles  darlingi  (1-7),  An.  pseudopunctipennis  (8),  and  An.  argyritarsis  (9)* 


Populations 

No.  of 

samples 

Mean  sample 
size/locus 

Mean  no.  of 
alleles/locus 

%  Polymwphic 
locit 

Mean  heterozygosity 

Direct-count 

HdyWbg  expected! 

1.  Belize 

4 

109.4 

2.0 

16.1 

0.063 

0.065 

(1.2) 

(0.2) 

(0.020) 

(0.020) 

2.  Amazonas,  Brazil 

4 

99.6 

2.5 

45.2 

0.101 

0.104 

(1.2) 

(0.2) 

(0.018) 

(0.019) 

3-  Sao  Paulo,  Brazil 

3 

93.4 

2.6 

48.4 

0.122 

0.130 

(1.2) 

(0.2) 

(0.023) 

(0.025) 

4,  Mato  Grosso,  Brazil 

1 

18.0 

1.9 

48.4 

0.122 

0.128 

(0.0) 

(0.2) 

(0.029) 

(0.029) 

5.  Venezuela 

3 

104.7 

2.7 

35.5 

0.098 

0.102 

(0.9) 

(0.2) 

(0.024) 

(0.024) 

6,  French-Guiana 

2 

93.5 

2.4 

19.4 

0.079 

0.084 

(0.6) 

(0.2) 

(0.021) 

(0.023) 

7.  Bolivia 

3 

23.2 

2.0 

45.2 

0.117 

0.125 

(0.5) 

(0.2) 

(0.025) 

(0.026) 

8.  An.  pseudopunctipennis 

1 

5.8 

1.3 

29.0 

0.081 

0.088 

(0.1) 

(0.1) 

(0.028) 

(0.029) 

9.  An.  argyritarsis 

1 

6.0 

1.2 

19.4 

0.048 

0.051 

(0.0) 

(0.1) 

(0.025) 

(0.021) 

*  Values  in  parentheses  are  standard  errors. 

t  A  locus  is  considered  polymorphic  if  the  frequency  of  the  most  common  allele  does  not  exceed  0.99. 
t  Unbiased  estimate.  HdyWjg  =  Hardy-Weinberg  equilibrium. 


(Table  4).  In  the  case  of  Idh-1,  which  has  four  alleles,  all 
populations  displayed  a  frequency  greater  than  72%  for  al¬ 
lele  100,  except  the  ones  from  Belize,  which  had  a  40% 
frequency  for  allele  100  and  59%  for  allele  86.  The  third 
locus  is  Pgm,  which  had  five  alleles.  The  Venezuelan  pop¬ 
ulations  had  a  frequency  of  70%  for  allele  80  and  28%  for 
allele  1(X),  whereas  all  the  other  populations  had  a  frequency 
greater  than  69%  for  allele  100. 

Regardless  of  the  allele  frequency  differences,  Nei’s  index 
of  genetic  identity  (above  the  diagonal)  indicated  a  high  de¬ 
gree  of  similarity  among  all  populations,  with  values  greater 
than  0.975  (Table  5).  The  phenogram  illustrated  that  An.  dar¬ 
lingi  populations  from  Brazil  and  Bolivia  had  high  genetic 
identities  (>  0.99).  Then,  by  decreasing  order  of  similarity, 
populations  from  French  Guiana,  Belize,  and  finally  Vene¬ 
zuela  followed  with  a  genetic  identity  of  0.981  (Figure  4A). 

The  comparison  of  An.  darlingi  populations  with  two  oth¬ 
er  anopheline  species  showed  that  with  An.  (Nys.)  argyri¬ 
tarsis,  the  genetic  distance  was  0.331,  and  with  An.  (An.) 
pseudopunctipennis,  the  genetic  identity  was  0.149  (Figure 
4B).  This  comparison  emphasizes  the  genetic  homogeneity- 
found  within  the  An.  darlingi  populations  studied.  The  pop¬ 
ulation  of  the  variant  paulistensis  presented  no  significant 
differences  compared  with  the  other  An.  darlingi  popula¬ 
tions. 

Random  amplified  polymorphic  DNA-poIymerase 
chain  reaction.  The  five  primers  used  in  this  study  produced 
35  bands  that  could  be  scored  in  An.  darlingi  and  17  in  An. 
nuneztovari,  including  four  shared  bands.  These  bands 
ranged  in  size  from  0.350  to  2.323  kilobasepairs  (Table  6). 
All  individuals  of  An.  darlingi  clustered  together,  as  did  all 
An.  nuneztovari,  each  with  100%  bootstrap  support  (Figure 
5).  Within  both  sjjecies,  there  was  some  evidence  of  geo¬ 
graphic  partitioning  of  most  populations.  Samples  of  An. 
darlingi  were  partitioned  into  four  main  clusters  (Figure  5). 
Cluster  I  comprises  most  (23  of  28)  of  the  Peixoto  de  Azev- 
edo  individuals  (central  Brazil);  Cluster  II  includes  samples 


from  Bolivia,  many  parts  of  Brazil,  French  Guiana,  and  Ven¬ 
ezuela;  Cluster  III  contains  many  (16  of  20)  of  the  Rio  Jacare 
Pepira  individuals  (southern  Brazil);  and  Cluster  IV  contains 
most  (9  of  10)  of  the  Belize  specimens.  Genetic  distances 
separating  all  clusters  of  An.  darlingi  were  small  in  com¬ 
parison  with  those  separating  the  two  clusters  of  An.  nunez¬ 
tovari,  which  included  samples  from  central  Brazil  (Peixoto 
de  Azevedo)  and  northeastern  Brazil  (Belem). 

Sequencing  of  1TS2.  Based  on  the  characterization  of  the 
rrS2  region  of  An.  darlingi,  no  variation  was  observed 
among  the  406  basepairs  of  the  ITS2  region  of  the  nine 
individuals  that  were  sequenced  in  the  present  study  from 
South  America,  including  Brazil  (Ituxi  River  and  Pereira 
Barreto)  and  French  Guiana  (Trou-Poisson  and  Maripasou- 
la).  However,  the  four  individuals  from  Belize  all  had  an 
identical  small  deletion  of  CCC  (Figure  6).  The  three  base- 
pair  differences  between  populations  from  Belize  and  South 
America  represents  a  value  of  0.74%  intraspecific  variation 
in  the  ITS2.  When  relationships  among  the  taxa  were  recon¬ 
structed  using  only  rrS2,  the  An.  darlingi  populations  from 
South  America  and  Belize  were  a  monophyletic  taxon  that 
was  not  allied  with  any  other  species  and  included  in  a  basal 
unresolved  polytomy. 

Morphology.  The  morphologic  study  showed  wide  vari¬ 
ation  in  the  pattern  of  pale  and  dark  spots  on  wing  costa, 
making  it  possible  to  differentiate  eight  classes  of  wing  pat¬ 
terns  (Figure  3A).  Costal  spots  are  characteristic  for  the  sub¬ 
genus  Nyssorhynchus,^  in  which  the  basal  dark  spot  (BD) 
is  absent  and  the  basal  pale  (BP)  and  prehumeral  pale  (PHP) 
are  fused  (BP  plus  PHP).  In  Type  I,  the  accessory  sector 
pale  (ASP)  is  present  and  divides  the  sector  dark  (SD)  into 
proximal  sector  dark  (PRSD)  and  distal  sector  dark  (DSD);’’ 
in  Type  11,  the  absence  of  the  ASP  resulted  in  the  fusion  of 
the  PRSD  and  DSD,  which  become  the  SD;  in  Type  HI,  the 
SP  and  ASP  are  absent;  in  Type  IV,  the  presector  pale  (PSP) 
is  absent;  in  Type  V,  the  subcostal  pale  (SCP)  is  absent;  thus, 
the  DSD  and  PD  are  fused;  in  Type  VI,  the  SP  is  absent;  in 
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Table  4 

Relative  allele  frequencies  for  three  loci  of  Anopheles  darlingi* 


Locust 

AJlele^ 

Belize 

Amazonas 

Brazil 

Sao  Paulo 

Mato  Grosso 

Venezuela 

French  Guiana 

Bolivia 

Got-l 

n 

114 

103 

99 

18 

107 

95 

24 

117 

— 

— 

— 

— 

0.009 

— 

— 

112 

— 

— 

0.01 

0.056 

0.023 

0.021 

— 

100 

0.982 

0.981 

0.616 

0.944 

0.808 

0.968 

0.958 

85 

0.018 

0.015 

0.374 

— 

0.056 

— 

0.042 

70 

— 

0.005 

_ 

— 

0.103 

0.011 

— 

(H) 

0.035 

0.039 

0.343 

0.111 

0.346 

0.063 

0.083 

Idh-l 

n 

102 

102 

99 

18 

107 

95 

24 

115 

0.02 

0.029 

0.056 

0.111 

0.084 

0.137 

0.042 

100 

0.392 

0.936 

0.889 

0.722 

0.752 

0.805 

0.896 

86 

0.588 

0.034 

0.051 

0.167 

0.164 

0.058 

0.063 

72 

— 

— 

0.005 

— 

— 

_ 

_ 

(H) 

0.52 

0.127 

0.222 

0.556 

0.439 

0.326 

0.208 

Pgm 

n 

111 

105 

101 

18 

108 

97 

24 

120 

— 

— 

0.005 

0.028 

— 

— 

0.021 

115 

— 

0.095 

0.089 

0.083 

0.009 

0.026 

0.167 

100 

0.797 

0.719 

0.851 

0.694 

0.282 

0.964 

0.729 

80 

0.203 

0.181 

0.05 

0.194 

0.704 

0.01 

0.083 

64 

— 

0.005 

0.005 

— 

0.005 

— 

— 

(H) 

0.315 

0.343 

0.277 

0.5 

0.426 

0.072 

0.458 

*  Dashes  represent  0;  bold  numbers  represent  the  allele  frequencies  discussed  in  the  Results, 
t  Got-J  =  glutamate  oxaloacetate  transaminase;  Idh-l  =  isocitrate  dehydrogenase:  Pgm  =  phosphoglucomutase. 
t  n  =  number  of  specimens;  h  =  heterozygosity  (direct-count)  per  locus. 


Type  VII,  the  PSP  and  SP  are  absent;  and  in  type  VIII,  the 
SP  and  SCP  are  absent. 

The  observed  frequency  distribution  of  wing  types  in  the 
different  populations  analyzed  varied  significantly  and  are 
population  dependent  (x^  =  161.34,  P  =  0.0016).  In  general, 
Type  I  was  the  most  frequent  type  of  wing  followed  by  Type 
VI.  In  the  Belize  population,  only  Type  I  and  n  were  present, 
while  Type  11  was  absent  from  the  populations  of  French 
Guiana,  Brazil,  and  Bolivia.  The  Venezuela  populations  of 
Caicara  and  Albarical  showed  up  to  six  different  t3fpes  of 
wings. 

Analysis  of  variance  (ANOVA)  of  the  proportions  of  the 
length  of  each  costal  wing  spot  divided  by  the  length  of  the 
wing  for  the  13  populations  showed  that  the  proportion  pre- 
humeral  dark/humeral  pale  (PHD/HP)  (Figure  3A)  is  a  valid 
diagnostic  characteristic  for  An.  darlingi  (F  =  78.1361,  P  = 
0.000001)  and  can  be  used  as  a  marker  to  differentiate  this 
species  from  other  morphologically  similar  ones  such  as  An. 
marajoara.  An.  braziliensis.  An.  argyritarsis,  and  An.  albi- 
tarsis. 

The  highest  Euclidean  distance  (ED)  was  observed  be¬ 


tween  the  Belize  and  Bohvia  populations  (0.14).  The  Boli¬ 
vian  population  is  separated  from  the  others  by  a  distance 
of  0.08.  The  Brazil  populations  of  Costa  Marques  and  Pe¬ 
reira  Barreto  are  similar  to  each  other  (ED  =  0.04)  and  to 
the  Venezuela  populations.  The  population  from  Colombia 
located  west  of  the  Andes  Mountains  is  very  close  (ED  = 
0.02)  to  the  population  of  San  Juan  de  Manapiare,  located 
in  the  Venezuelan  Amazon,  east  of  the  Andes  (Figure  2). 

Analysis  of  the  morphology  of  the  proboscis  and  legs  re¬ 
sulted  in  vahd  diagnostic  characteristics  that  could  separate 
A«.  darlingi  populations  such  as  proboscis  (P)  and  forefemur 
(P/Fe-I),  pale  spot  on  foretarsomere  (PS-I,/Ta-I,,  PS-l2/Ta-l2, 
PS-l3/Ta-l3)  and  on  hindtarsomere  (PS -1112714-1112)  (Figure 
3B).  Galvao  and  others"  pointed  out  that  the  variant  paulis- 
tensis  had  more  dark  scales  on  the  hindtarsomere  2  than  the 
An.  darlingi  described  by  Root.  Nevertheless,  for  this  char¬ 
acteristic,  we  found  that  the  Pereira  Barreto  population  (pau- 
listensis  type  locality)  was  similar  to  the  Bolivian  popula¬ 
tion.  The  ANOVA  of  these  proportions  for  An.  darlingi.  An. 
marajoara.  An.  braziliensis.  An.  argyritarsis,  and  An.  albi- 
tarsis  was  significantly  different  {P  <  0.(KX)1),  and  the  best 


Table  5 

Matrix  of  genetic  distance  of  Anopheles  darlingi  (1-7),  An.  pseudopunctipennis  (8)  and  An.  argyritarsis  (9)* 


Population 

2 

3 

4 

5 

6 

7 

8 

9 

1.  Belize 

8c  9k  8c 

0.987 

0.978 

0.988 

0.983 

0.985 

0.980 

0.114 

0.307 

2.  Amazonas,  Brazil 

0.014 

**** 

0.992 

0.995 

0.987 

0.992 

0.993 

0.134 

0.334 

3.  Sao  Paulo,  Brazil 

0.022 

0.008 

0.990 

0.976 

0.984 

0.992 

0.131 

0.341 

4.  Mato  Grosso,  Brazil 

0.012 

0.005 

0.010 

sksksfcsk 

0.986 

0.988 

0.991 

0.125 

0.337 

5.  Venezuela 

0.017 

0.013 

0.024 

0.014 

8c  8c  8c  8c 

0.979 

0.979 

0.124 

0.322 

6.  French  Guiana 

0.015 

0.008 

0.016 

0.012 

0.022 

8c  8=  8c  8c 

0.985 

0.126 

0.313 

7.  Bolivia 

0.020 

0.007 

0.008 

0.009 

0.021 

0.015 

8: 8: 8c  8c 

0.118 

0.367 

8.  An.  aygyritarsis 

1.181 

1.098 

1.076 

1.087 

1.134 

1.161 

1.003 

8c  8=  8s  8c 

0.324 

9.  An.  pseudopunctipennis 

2.172 

2.007 

2.032 

2.081 

2.089 

2.069 

2.140 

1.128 

8: 8c  8c  8c 

*  Values  above  the  diagonal  are  Nei’s  unbiased'^  genetic  identities;  values  below  the  diagonal  are  Nei’s  unbiased"  genetic  distances. 
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Figure  4.  Unweighted  phenogram  from  Nei’s‘“  unbiased  genetic  distance  matrix.  A,  for  all  Anopheles  darlingi  populations  (cophenetic 
correlation  =  0.815).  B,  for  all  An.  darlingi  populations.  An.  pseudopunctipennis,  and  An.  argyritarsis  (cophenetic  correlation  =  0.983). 


Table  6 

Approximate  sizes  of  random  amplified  polymorphic  DNA  (RAPD) 
fragments  used  in  comparison  of  Anopheles  darlingi  populations 
and  comparison  of  An.  darlingi  with  An.  nuneztovari 


Primer  Size  of  RAPD  fragments  (kilobasepairs)* 


B12  1.569,  1.501t,  1.226,  1.140,  1.135,  0.940,  0.911,  0.892t, 

0.888,  0.861t,  0.705t  0.675,  0.552,  0.439*,  0.423t 
A09  1.888,  1.593,  1.323**,  1.071,  0.872,  0.817*,  0.781, 

0.449,  0.432*,  0.350 

B04  1.150*,  0.760,  0.707,  0.468*,  0.429,  0.422*,  0.409,  0.358 

A07  2.322*,  1.721*,  1.489,  1.297,  1.152,  1.045*,  0.976*, 

0.790,  0.735,  0.636,  0.589,  0.509 
CIS  1.935*,  1.252,  0.718 

Occurrence  of  fragments:  no  annotation  =  An.  darlingi. 

+  An.  nuneztovari. 

$  Both  species. 


diagnostic  characteristic  to  separate  all  four  species  is  the 
proportion  PS-in2/Ta-IIl2.  This  characteristic  was  previously 
used  by  Linthicum^  to  separate  An.  albitarsis  and  An.  ma- 
rajoara  only. 

DISCUSSION 

Based  on  the  isozyme  analysis,  the  mean  ±  SE  hetero¬ 
zygosity  of  0.100  ±  0.023  for  An.  darlingi  (index  that  re¬ 
flects  the  genetic  variabUity)  was  similar  to  the  average  value 
found  in  other  Diptera  (0.115  ±  0.(X)9).“  This  species 
showed  a  higher  genetic  variability  than  the  An.  pseudo¬ 
punctipennis  populations  (0.059  ±  0.020)  studied  by  Man- 
guin  and  others.*^  However,  lower  genetic  variability  was 
found  for  the  BeUze  populations  (0.063)  compared  with 
South  American  populations  (0.079-0.122,  with  a  mean  ± 
SE  of  0.107  ±  0.017).  Belize  and  other  countries  of  Central 
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Figure  5.  Phenogram  showing  relationships  among  124  individuals  of  Anopheles  darlingi  and  13  An.  mmeztovari.  Branch  lengths  are 
proportional  to  1  -  S,  and  the  scale  bar  is  equal  to  0.1.  On  selected  branches,  the  bootstrap  value  of  100  repetitions  is  shown.  The  roman 
numerals  indicate  the  four  primary  clusters. 


America  are  part  of  a  rather  narrow  and  restricted  land  mass 
that  is  scored  by  many  short,  parallel  river  systems.  As  a 
consequence,  populations  of  An.  darlingi  from  Belize,  and 
more  globally  from  Central  America,  have  a  restricted  dis¬ 
tribution  compared  with  the  large  land  mass  and  extensive 
river  systems  occurring  in  South  America.  It  may  be  that  the 
loss  of  population  heterozygosity  in  Belize  reflects  the  nor¬ 
mally  expected  loss  of  polymorphic  alleles  in  populations 
that  are  confined  to  restricted  and  isolated  environments.  In 
other  words,  the  gene  flow  of  An.  darlingi  populations  in 


Central  America  occurs  in  a  more  reduced  area  than  in  South 
America,  causing  a  bottleneck  effect  that  resulted  in  less 
polymorphic  populations.  Our  hypothesis  is  that  a  northern 
extension  of  populations  from  South  America  established 
An.  darlingi  populations  in  Central  America.  It  may  also  be 
that  the  modem  discontinuous  geographic  distribution  of  this 
species  occurred  secondarily  to  the  original  northward  ex¬ 
tension  of  An.  darlingi  populations. 

Overall  genetic  identity  within  An.  darlingi  populations 
was  high  (0.976-0.995),  and  these  identities  are  markedly 
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Figure  6.  Sequence  of  34  nucleotides  of  the  5.8S  region,  the  entire  406  nucleotides  of  the  internal  transcribed  spacer  2  (ITS2)  region,  and 
24  nucleotides  of  the  28S  region  for  Anopheles  darlingi.  The  arrowheads  represent  the  first  and  last  sequences  of  the  ITS2  region.  SA  = 
South  America;  BE  =  Belize,  which  differs  only  by  a  small  three-base  deletion  (CCC)  at  position  179-181.  The  sequence  has  been  deposited 
in  GenBank  (accession  number  U92337). 


higher  than  the  value  of  0.85,  suggested  by  Avise*’  as  the 
lower  limit  for  conspecific  populations.  For  instance,  in  other 
anopheline  species,  values  ranged  from  0.84  to  0.88  between 
both  species.  An.  gambiae  s.s.  and  An  arabiensis,’’^''^  and 
from  0.92  to  1.0  among  An.  pseudopunctipennis  popula¬ 
tions.''^  In  the  study  of  Steiner  and  others,^®  an  An.  darlingi 
population  from  southern  Brazil  had  values  ranging  between 
0.592  and  0.772,  suggesting  the  existence  of  more  than  one 
species.  In  comparison  with  another  study,'*  genetic  identity 
of  southern  and  northern  Brazil  populations  has  been  found 
to  be  much  higher,  with  a  value  of  0.952  within  the  range 
of  intrapopulational  variation.  Similarly,  the  range  of  our  ge¬ 
netic  identity  values  for  all  the  Brazil  populations  was  both 
narrower  and  higher  (0.976-0.992)  than  the  former  one. 

Although  some  moderate  differentiation  was  found  in 
three  loci  {Got-1,  ldh-1,  and  Pgrri),  the  mean  Fst  values 
showed  negligible  differentiation  among  An.  darlingi  and  no 
fixed  differences  were  found.  In  addition,  the  moderate  mean 
F[s,  with  a  value  of  0.039,  suggests  that  random  mating 
among  the  populations  of  An.  darlingi  is  occurring.®^  Com¬ 
parison  of  the  isozyme  analysis  between  An.  darlingi  and 
An.  pseudopunctipennis*^  showed  a  higher  homogeneity 
among  populations  of  the  former  species,  with  Fst  values  of 
0.102  and  0.375,  and  Nei  genetic  identities  ranging  between 
0.976  and  0.995,  and  0.924  and  1.0,  respectively.  In  the  case 
of  An.  darlingi,  there  were  no  clear  groupings  of  populations 
as  found  for  An.  pseudopunctipennis. 

Our  RAPD-PCR  results  support  the  conclusion  that  the 
widespread  sampled  populations  of  An.  darlingi  are  conspe¬ 
cific.  In  addition  to  high  similarity  values  among  the  popu¬ 
lations,  comparisons  with  An.  nuneztovari  demonstrated  the 
level  of  differentiation  expected  between  species.  Of  the  35 
markers  used  in  this  analysis,  seven  are  monomorphic,  i.e., 
occur  in  >  95%  of  the  individuals  tested,  and  two  of  these 
occur  in  all  individuals.  Wilkerson  and  others^’^^  found  sim¬ 


ilar  numbers  of  correlated  diagnostic  markers  in  each  of  the 
four  species  of  the  Albitarsis  complex  (3-7)  using  the  same 
60  RAPD  primers  as  in  the  present  study. 

Our  analysis  produced  four  large  and  three  small  clusters, 
indicative  of  geographic  partitioning.  Since  An.  darlingi 
clusters  usually  include  individuals  from  widely  separated 
populations,  e.g.,  cluster  III  includes  individuals  from  central 
and  northern  Brazil  and  Bolivia,  the  case  can  also  be  made 
that  the  level  of  mixis  implied  by  the  clustering  would  likely 
be  sufficient  to  ensure  panmixis  (random  mating). 

Interestingly,  the  two  clusters  in  the  RAPD-PCR  data  for 
An.  nuneztovari  (Bel6m,  northeastern  Brazil  and  Peixoto, 
central  Brazil)  are  similar  to  the  distinctive  mitochondrial 
DNA  lineages  of  12  An.  nuneztovari  populations  in  Bel6m 
and  central  Amazonian  Brazil  (Manaus  and  Boa  Vista),®*  un¬ 
derscoring  some  congruence,  in  this  case,  between  two  dif¬ 
ferent  molecular  markers. 

Intraspecific  variation  in  the  ITS2  region  of  21  members 
of  the  Nyssorhynchus  subgenus  varies  from  0%  to  2.8%. 
Therefore,  the  value  of  0.74%  for  An.  darlingi  represents  a 
low  value  for  intraspecific  variation  of  the  subgenus.  Most 
of  this  variation  is  in  the  form  of  insertion-deletions  involv¬ 
ing  a  few  base  pairs  (as  is  the  case  for  An.  darlingi  from 
Belize  in  the  present  smdy)  or  single  base  point  mutations. 
The  identical  ITS  sequences  across  all  of  South  America 
suggest  a  single  panmictic  breeding  population. 

We  conclude  that  samples  from  Mexico,  as  well  as  other 
Central  American  populations  and  those  from  Colombia, 
would  have  to  be  sequenced  to  determine  whether  the  small 
deletion  in  the  Belize  population  represents  isolation  by  dis¬ 
tance,  immigration,  or  extinction  events.  Isolation  by  dis¬ 
tance  could  result  from  Belize  populations  being  near  the 
northern  end  of  a  north-south  cline  or  could  be  due  to  im¬ 
migration,  or  extinction  events  could  be  related  to  the  re¬ 
ported  discontinuity  of  An.  darlingi  populations  in  Central 
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America.  Both  distance  explanations  seem  to  be  relatively 
weak  since  there  is  no  evidence  of  clinal  variation  (for  ex¬ 
ample,  the  rrS2  sequence  is  identical  between  El  Juval  in 
western  Venezuela  and  Sao  Paulo).  If  the  Belize  populations 
had  been  independently  evolving  for  a  long  period  of  time 
due  to  an  historical  vicariant  (of  different  origin)  event,  var¬ 
iance  in  the  rrS2  from  Belize  compared  with  populations  in 
South  America  should  be  greater.  Therefore,  we  suggest  that 
the  small  deletion  in  the  Belize  populations  of  An.  darlingi 
is  most  likely  the  result  of  a  recent  introduction  event  caused 
either  by  human  action  or  some  other  extrinsic  force.  Future 
studies  could  test  this  hypothesis  by  determining  which 
South  American  populations  have  the  greatest  similarity  with 
the  Belize  populations  and  whether  other  Central  American 
and  Mexican  samples  are  more  similar  to  populations  from 
Belize  or  from  South  America. 

The  study  of  head,  legs,  and  wings  of  17  populations  of 
An.  darlingi  showed  wide  variation  in  morphologic  charac¬ 
teristics.  On  the  wing  costa,  eight  different  patterns  were 
described  based  on  the  presence  or  absence  of  one  or  more 
spots.  Geographic  variations  in  the  pattern  of  the  wing  costa 
have  been  reported  for  other  Nyssorhynchus  such  as  An.  al- 
bitarsisA’  and  An.  nuneztovari.^^  The  Euclidean  distances 
separating  groupings  of  diagnostic  wing  characters  suggest 
that  gene  flow  occurs  among  populations  across  all  geo¬ 
graphic  locations  of  populations  included  in  these  analyses. 
A  previous  study*’  on  Brazil  populations  showed  that  the 
variations  were  intraspecific  “because  no  patterns  were  pop¬ 
ulation  specific”. 

Based  on  isozyme,  RAPD-PCR,  ITS2  sequence,  and  mor¬ 
phologic  analyses,  all  An.  darlingi  populations  collected 
from  most  parts  of  the  species’  range  in  Central  and  South 
America  demonstrated  a  genetic  similarity  that  is  consistent 
with  the  existence  of  a  single  species.  While  some  genetic 
and  morphologic  variations  were  found  among  the  popula¬ 
tions  surveyed,  no  complete  barriers  to  gene  flow  were  iden¬ 
tified.  The  relative  homogeneity  of  An.  darlingi  populations 
might  indicate  that  the  isolated  population  of  Central  Amer¬ 
ica  could  have  originated  from  South  America.  As  presently 
defined,  the  limited  range  of  An.  darlingi  in  Central  America 
suggests  a  recent  introduction  of  the  species.  This  limited 
distribution  also  raises  the  possibility  that  the  species  could, 
over  time,  expand  into  or  reinvade  more  Central  American 
countries  and  more  states  in  southern  Mexico.  The  former 
hypothesis  could  be  tested  using  microsatellites,  RAPD- 
PCR,  or  ITS2  sequences  for  the  characterization  of  An.  dar¬ 
lingi  populations  throughout  southern  Mexico,  Central  and 
South  America  to  determine  more  precisely  the  phylogeny 
of  the  species. 
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